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A theoretical approach to the influence of one-dimensional lattice fluctuations on electronic 
properties in weakly localized spin-Peierls systems is proposed using the renormalization group and 
the functional integral techniques. The interplay between the renormalization group flow of corre- 
lated electrons and one-dimensional lattice fluctuations is taken into account by the one-dimensional 
functional integral method in the adiabatic limit. Calculations of spin-Peierls precursor effects on 
response functions are carried out explicitely and the prediction for the temperature dependent 
magnetic susceptibility and nuclear relaxation is compared with available experimental data for 
(TMTTF)2PF6. 



I. INTRODUCTION 



OO ' Among the possible structural phase transitions in low-dimensional systems, the spin-Peierls (SP) instability has 
recently received particular attention. [Q,^. In contrast with the Peierls transition, whose normal phase is metalhc, 
pronounced insulating properties and antiferromagnetic spin correlations are found far ahead of the SP ordering 
I temperature. However, owing to the quasi-one-dimensional character of these systems, precursor effects should appear 
in both types of transitions. Since spin correlations are essentially one-dimensional in character, lattice fluctuations 
■ should emerge as diffuse scattering in x-ray experiments, and as a so-called pseudo-gap effect in magnetic properties. 
' The pseudo-gap effect is well known to be the hallmark of strongly anisotropic Peierls systems , but surprisingly, 
its observation is a rather uncommon feature of most SP systems and the underlying explanation for this is not really 
known so far ||^. Correspondingly, the incitement to achieve a theoretical description of fluctuation effects in the 
SP case did not receive any experimental justiflcation, so that most theoretical efforts were rather devoted to the 
d ' description of the ordered state However, recent observations of SP fluctuation effects in the organic systems 

g ! (TMTTF)2X and (BCPTTF)2X with X^PFg, AsFg, has prompted a renewed interest for this problem In 
both cases, one-dimensional SP precursors are clearly seen in structural and magnetic measurements well above the 
true SP ordering temperature. 

This paper will be entirely devoted to the theoretical description of one-dimensional lattice and electronic corre- 
Q lations in SP systems of the (TMTTF)2X series. Unlike (BCPTTF)2X, which are already good insulators at room 
temperature, (TMTTF)2X are weakly localized and still show a metallic behaviour in this temperature range. It 
is only around Tp ~ 220K that carriers undergo a Mott-Hubbard localization, below which the resistivity becomes 
thermally activated ||l^. Albeit the magnetic susceptibility remains unaffected at Tp, NMR spin- lattice relaxation rate 
measurements have revealed the existence of strong ID antiferromagnetic spin correlations below this temperature 
range [tHi^. Below 60K, one-dimensional SP lattice softening is horned out by x-rays and has a marked influence 
on the temperature dependence of both susceptibility and nuclear relaxation . In order to tackle the problem of 
one-dimensional electronic and lattice fluctuations effects in SP systems we shall combine the renormalization group 
technique for fermions and the functional-integral method ||9|,p^ . A similar approach has been recently applied -with 
some success- to strongly localized SP systems of the (BCPTTF)2X series [Q. Since the itineracy of carriers is much 
more pronounced in the high temperature domain of (TMTTF)2X series, both spin and charge degrees of freedom 
must be considered in the renormalization flow, which in practice makes the extension of the previous approach 
slightly more complicated. 



II. THE MODEL AND THE PARTITION FUNCTION 



In setting up the essential ingredients of the model, it is first convenient to study the partition function of the 
electron-acoustic-phonon interaction model, namely when direct electron-electron interaction is absent. The corre- 
sponding hamiltonian is well known to take the form p3[ : 
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p,k,a Q {p,k,Q,a} 

where ep(fc) — vp{pk — kp) is the hnearized spectrum around the ID Fermi points pkp = ikp for free fermions and 
'^(Q) = sin((5a/2)l is the dispersion relation for the phonons (h = 1). Here vp is the Fermi velocity and lou is 
the Debye frequency, which coincides with the phonon energy at Q = ±.2kp{— ztn/a, when the band is half-filled). 
As for the interacting part, the electron- phonon coupling constant is given by: 

5(fc, Q) = z4-=^== sin(Qa/2) cos(fca + Qa/2), (2) 

where M is the molecular mass and ga — —dta/dx > is the spatial modulation of the longitudinal hopping integral 
ta- The dynamics of phonons being purely harmonic, their integration in the partition function Z = Tve~^^ can be 
carried out exactly, yielding an effective retarded interaction among electrons (/3 = 1/T, fcs = 1). The result can be 
put in the following functional-integral form 

Z= j [dr][d^] e^°['^*''^l exp(- dTidr2 W(ri - T2)), (3) 

over the anticommuting Grassman variables ■!/)'•*-'. The free electronic action is given by 

S°[r,^] = f dT E'A;..(A;,r)[GO(fc,T)]-Vp,<.(fc,T), (4) 
JO , 

where Gp{k,T) ~ [—d/dr — ep{k)]~^ is the free electron propagator. As for the retarded interaction W(ri — T2), we 
will proceed to a "g-ology" decomposition which only retains 2kp phonon exchange between electrons at ikp, with 
the result 

H(ri-T2)= (2L)-i J2 9i,ph{Ti~T2) i;;^,^{ki+p2kp + q,Ti)^*_p^^^{k2-p2kp-q,T2) 

{p,k,q,a} 

X '0p,cr2(fc2,T2)'0_p,ai(^l,n) 

+ (2^)"^ E 53,ph(Ti - T2) Vlp,ai(^i +p2kp - q - pG,Ti)iptp,^^{k2 -p2kp + q,T2) 

{p,k,q,a} 

X '(/'p,<T2(fc2,T2)i/'p,<Ti(A:i,ri), (5) 
where G — Akp is a reciprocal lattice vector. The retarded backscattering and umklapp coupling constants are given 

by 

5i,ph(Ti - T2) = g{-pkp ,p2kp)g{pkp , ~p2kp)D° {2kp,Ti - T2) 
- ^^'^ D°i2kp,n-T2), 



MujD 

g3,ph(ri - T2) = -?y5i,ph(ri - T2), 

where D°{2kp,Ti - T2) = -[e-'^o\ri-T2\ _^ 2{e'^'^° - l)-icosh(wD(Ti - T2))] is the phonon propagator at 2kp. Here, 
ry < 1 is a positive constant that reflects the half-filled character of the band. For a single half-filled band, 77 = 1 while 
for systems like (TMTTF)2X with a quarter-filled band and a small dimerization gap, one can take 77 ^ 1 From 
above, one observes that, for exchange of 2kp acoustic phonons, 53. ph and gi.ph are opposite in sign. The retarded 
interaction can then be written uniquely in terms of composite fields: 

- T2) = / dx E.9ph(n - T2)0''*{X,T,)0'\X,T2), (6) 

where we have introduced the composite field 

0^^*(x,T) = i(0*(x,r)+MO(x,T)), 

with 

0{x,t) = ^1p*_,^{x,T)lp+^a{x,T). 
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These corresponds to "site" (M = +) and "bond" [M = — ) charge-density-wave (CDW) correlations at half-fiUing. 
The related combinations of couplings are given by 

Qpli'Tl - T2) = ,gi,ph(Tl - T2) + A/53,ph(Tl - T2). 

In order to analyze these correlations, we apply an Hubbard-Stratonovich transformation on the retarded part of the 
interaction 

Z = I [d(^*]m e'"^^*'^^ exp{- / dxdr^dr^ ^ 0^^(:r,Ti) | 5ph(n - T2) cj,"{x,T2)} 

M=± 

xexp|- / dxdr AAfO^'^la;, t)0^^(x, t)}, (7) 



M=± 

where 0^ are real auxiliary fields for site and bond fluctuations respectively, and A+ = 2, A_ = 2i. It is worth noting 
that, in the absence of umklapp scattering, g^^^ = g^y^ and there is no difference between bond and site density wave 
fluctuations so that both auxiliary fields will combine to produce amplitude and phase fluctuations of the complex 
field (j) = <jy^ -\- i(f)~ =101 e^'^ , as found in incommensurate Peierls systems. 

Unretarded repulsive electron-electron interaction at half-filling will promote bond with respect to site electronic 
correlations. These essential ingredients of the SP instability can be added to the model by following the "g-ology" 
prescription p5[ , in which the direct interaction among right- and left-moving carriers is decomposed in terms of 
backward (gi), forward (172) and (33) umklapp scattering coupling constants. The full partition function of the model 
in the Fourier-Matsubara space then becomes 

z[h^'] = J [dr][dndci^''] exp{s°^[r,^] + s°[ci>^'] + sx[r,i^,cf>^'] + si[r,^] + Sh[r,i',h'']} 
[dr][dndci^''] exp{ [G°(fc)]-v;,.(fc)V'i.,a(fc) - E I '('^'(Q) I'l 9''m r' 

p,k,a Q 

Q 

+ 2Z ^ 32-0^,^1 (fcl + q)'^%.a2 (^2 - q)^p,a^ {ki)'lp-p,a2 (h) 

2T ^ ~ ~ 

+ — 2^ gs't/jp.aAki + Qp)'0p,<T2 (^2 - Qp +pG)?A_p^ff,(A:i)V'-p.<T2(^2) 

{p,k.Q,a} 

+ E(of(^)'^rw)+H-c)}' 

where k = (k,Lj„ = {2n + l)TrT)), Qp = {pQ,uj„, = 2^mT), G = (4fcF,0), C^pik) = [iuj^ - ep(k)]-\ and g^\u„,) = 
8gl{MujDy\l - M-q)D°{2kF,uJm)- The integration measures are [d'ip*][d'ip] = U.p%„dip* ,^{k)d'ipp.„{k) and [d0^^] ^ 
^MQ>o I dphi^-in) 1^^ d(j)^'^ {Q)d(j)^^ {—Q) for fermion and auxiliary fields, respectively. Here 

D'^{2kF,u;m)^^^^^, (9) 

corresponds to the bare phonon propagator. Finally, Sh is an additional term which couples an infinitesimal source field 
h^^ to site bond (Af = ±, ^ = 0) CDW and site (M = +, ^ ~ 1,2,3) SDW correlations. These correspond to composite 

fields defined by 0^'*iQ) = i(0;(Q) + MO^(O)), with 0;.o,i.2,3(Q) = iT/L)^/^Ez,^,pr+Jk + Q)af^.,pik), 
"■1,2,3 and tJo being the Pauli and the identity matrices respectively. In the following the source-field term will be 
useful for the calculation of relevant response functions for the SP instability. 
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III. RENORMALIZATION GROUP RESULTS 



The problem of low-energy electronic and lattice correlations reduces to the study of interacting electrons coupled 
to the fluctuating field 0*^, which can be done by first applying the renormalization group approach developped in 
Ref. [|l6l,f0| for the fermion degrees of freedom. It consists in first integrating high-energy fermion states, namely we 
write V^^—> '0''*^ + where 'i/S^*-' describes degrees of freedom to be integrated over in the outer energy shell of 
tickness ^Eo{£)d£ on both sides of the Fermi level at ikp and for all cUn- Here Eo{£) = Eoe~^ is the band energy 
cut-off at the step £ and Eq = 2Ef is the initial band width, which is twice the Fermi energy. Keeping the (/)^'^'s fixed, 
this is formally written as 



n 



< 

which leads to a recursion relation for electronic parameters of the action including those related to source fields 
for the calculation of response function (see below) and for T[4>^'^]^, which is a free-energy functional which collects 
all the contributions in the auxilliary field (j>^'^ . At finite temperature, the partial integration is conducted down to 
£t = ln(i?F/2^)- In the present renormalization group procedure, electronic degrees of freedom are treated in the 
continuum limit so that the discreteness of the underlying lattice will be neglected in the following . 

A. Electronic part 

1. Couplings and one-particle self-energy 

It is useful for the following discussion to recall the well-known two-loop RG results for the purely electronic part 
(0*^ 0) in zero field {h^ = 0). The effect of the lattice fluctuations and the calculation of relevant response func- 
tions will be considered afterwards. In this purely electronic limit, the single-particle fermion propagator transforms 
according to the recursion relation zi{£ -\- d£)[Gp]~^ = zi{£)zi{d£)[Gp]~^ . At two-loop level the outer-shell correction 
to Zi{d£) comes from the n = 2 (Sj )o.s terms of (p^, which yields the flow equation: 



dln^f^ 1 



:{[252W-giWp-5|W+3[5?-h35|]}, (11) 



d£ 16 

where the gi{£) = gi{£){'KVp)~^ are the normalized coupling constants at I. The latter transform as 'gi(£ + d£) = 
gi{£)zi'^ {d£)z2,3,4{d£) , which, together with zi, are obtained from one-loop (Sf)o.s and two-loop (<5f )o.s outer shell 
corrections to two-particle four-point vertex functions ri^2,3(^ + d£) — Z2.3^4{d£)Ti_2,3{£)- This is known to yield the 
flow equations 0,|l|,|l§ 

dgi_ ^2 1-3 
'd£~ 2^1' 

d{2g2-gi) ~2m 1/o~ ~m 
^= 33(1 --(232-51)], 

^= 53(252 - gi)[l - \{2g2 - 51)] - \9l (12) 

For (TMTTF)2X compounds, 51 ~ 52 > are repulsive at ^ = 0, whereas a finite dimerization of the TMTTF stacks 
can be parametrized by a small and positive 53 <^ 51 |lj,^ . One then gets the inequality gi — 2g2 < 53 indicating that 
umklapp scattering is relevant. The flow of charge couplings 252 — gi and 53 then scales towards strong coupling while 
the spin coupling gi is found to be marginally irrelevant with the fixed point values g^{£ — > 00) — > 2, 52 ~* ^) ~* 
and gl — + 0. The strong coupling sector of the RG flow is reached when 53 (£ = £p) « 1, where £p — \n{Ep/Tp) defines 
the temperature scale for the presence of a charge (Mott-Hubbard) gap Ap = 2Tp. Neglecting transients between 
weak and strong coupling regimes, one-particle self-energy corrections resulting from the solution of Eq. ( pi] ) can be 
put in the following scaling form 
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^^\^)«%-\^.)(£^oW/Apr^ (13) 

which depicts the reduction of the quasi-particle weight at the Fermi level with 0* (91,92,93) = 3/4 being evaluated 
at the fixed point. 

2. Response functions 

The pertinent electronic response functions involved in the description of the SP instability in one dimension are 
those related to site SDW (/x ^ 0,M = +) and bond CDW (/x = 0, M = -) correlations Q. They can be readily 
calculated via the renormalization of the source field term 5*^, which, at scale £, reads 

S,[r,i^,h'%^ {z^\e)hf,'%Q)0'J{Q)+U.c.-xf:'{2kF,£)hf!*\Q)hf,'{Q)}, (14) 

where zjf {£) is the renormalization factor for the pair vertex part, while 

Xf (2fcF,£) = -(nvF)-' fxf:'inde', (15) 
Jo 

is the 2kp response function and xff = [zffY is the auxiliary susceptibility in the (/i, M) channel considered . The 
one- and two-loop corrections to z^^ come from the outer shell averages {ShSi)o.s and {ShSj)o.s for n — 2 and n = 3 
respectively. This leads to the flow equation 

dlnv*^ 1 1 

= g^'W " 2 {5?(^) + glii) ~ ~9im2{£) + i^rS)), (16) 

where g±(^) ^92{t}T93.{£)-2gi{t} for the CDW channel and 5+_^o (^) =52(^)+53(^) for the SDW channel. Following 
the example of (^^ for the one-particle self-energy renormalization, can be expressed in the following scaling form 

xf(^)«xf(^p)(i?o(^)/Ap)-^- , (17) 

where x*^(^p) is the weak coupling contribution below the charge gap. The exponents 7*^3 = 7*^ = 3/2 obtained 
by evaluating the r.h.s of ( [l^ ) at the fixed point indicate that only site SDW and bond CDW correlations are singular 
at half-filling. It should be mentionned, however, the g* obtained at two-loop level are well known to overestimate 
the values of the exponents. Higher order contributions are expected to bring them closer to the exact value 

711,3 = 7*" ^ 7* - 1, (18) 

which is known to result from more elaborate calculations at ^ ^ [|l^ . 

Before closing this section, it is useful for applications (see §111. B. 5 and §IV.B) to compute the imaginary part of 
the retarded response function at low (real) frequency. From Ref. [ p^ , it is directly connected to as follows: 

Imxf(2fcF + g,a;)=xf(7^)Inix"('? + 2fcj.,c^), u ^ Q (19) 

for vpq ^ r, where 

2T 

Imx''(g + 2fcF,c^) = — Im^ G'i(fc,w„)G°(fc + 2fcF + g,cc;„+cj + iO+) 

= TT. W^O, 20) 

is the imaginary part of the low-frequency non-interacting response function near 2kp (see Appendix A. 2). 
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B. Influence of lattice fluctuations 



1. Ginzburg- Landau functional 

When the lattice auxihary field is taken into account, the partial integration (||) generates a series of terms in 
n > 2 powers of the 0^^'s via the closed fermion loops {S^)/nl. This will not only lead to corrections for S^[(j)^^]i 
for n — 2, but it yields a recursion relation for the quantum Ginzburg-Landau free-energy functional ^-"[0*^]^, when 
combined to the n > 2 terms [^^. Since only 2fci? bond correlations are singular in the CDW channel for repulsive 
couplings, only the dependence on the bond field 0^ needs to be retained, with the result 

/3^[<^-],+<i, = /3^[0-], - ^^B„(d^) 0-(gi)...0^(g„)jj._g. (21) 

n>2 Q 

The quadratic {n — 2) and quartic {n — 4) outer shell contributions lead to 

B2{d£) = (7rup)~i[z-(^)]2 d£, 

. (22) 

where C(3) ~ 1.2 . . .. For the problem at hand, one can assume that adiabaticity between electrons and the lattice 
is sufficiently strong that quantum lattice corrections can be neglected, thereby allowing the static limit. Up to the 
quartic term, one obtains the following Ginzburg-Landau free energy functional (in real space) 

T[r]tr= j dx[a{T){<p-{x)f +c(^^ +6(0-(x))'], (23) 

at temperature T, where a(T), c and b are the Ginzburg-Landau parameters. The coefficients a(T) and c of the 
quadratic term are respectively given by 

«(r) = [5ph(0)]-'+X"(2fcF,T) 

«a'(r/r^p-i) 

c « ^a'ivF/TTTlpf, (24) 



where a' = {TrvF)~^X~{Tsp)- Here the bond CDW response function x {'^kp^T) in Eq. (p^, together with (p^), 
have been used in the linearization of a(T) around the SP mean-field temperature 

Tlp^T,{\g;^\x-iT,))'^''- (25) 

for Tp 3> Tgp. It is worth noting that since | g~j^ \ cc g'^, then Tgp oc g^, which turns out to be the same power 
dependence on the electron-phonon coupling constant that was found in previous mean-field calculations on more 
localized systems As for the rigidity parameter c, it is readily obtained in Fourier space from the Q = 2kF + q 
expansion of X~ (Qi ^sp) using the approximate expression 



e{Tlp, 2kp + 9) = In |^[1 + {vFq/nTlprr^^ , 



as boundary condition of the logarithmic integration in (p^ ) near Tgp. This form essentially coincides with the Q 
dependence of the elementary Peierls bubble near 2fcp. Finally, the flow of the mode-mode coupling in (^2|) is stopped 
at ix" and reads 



SP 



for Tp > TOp. 



Fluctuation effects in (f)" below Tgp are then governed by the classical functional integral Z — 
J[d(j)~'] exp{— (3J-[(j)~]) , which can be carried out exactly using the transfer matrix method pl[ . An important quantity 
to compute is the static ID SP response function xsp- From the results of Ref. p^, one finds 
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0.0 0.4 0.8 





T/Tsp 

FIG. 1. Temperature variation of the mean square of the SP order parameter (left scale ) and the inverse of the correlation 
length (right scale) normalized by the lattice constant. 



_ 2/3 (( (r )^ )) gsp 



1 + 

where 



(( (...) )) = z-i [dr] {...)eM-mr]), 



denotes a statistical average over ip . Here ^sp is the correlation length of the real order field (p , which grows 
exponentially below Tgp according to 



where 



^SP^-.f^e"/^, (28) 
p\ \ a\ 



P=^^^^^-3.8T^„ (29) 



for T <C Tgp. The temperature profile of both £^sp and the mean square fluctuations {{{(f> )^)) are plotted in Figure 
1. 



2. Pseudo-gap effect 

The infiuence of the static SP lattice fluctuations governed by (|2^) on one-particle electronic properties at the step 
£t of the RG procedure will be analyzed by considering to leading order the one-electron self-energy contribution of 
Figure 2 jsj, which can be obtained from (^) after an integration over the 4>~ field. The result is 

W-}) = -TL-h^'z^ J2 - P^'^F - 9. *^n)xsp(2A:F + q, T). (30) 

q 

The bare propagator appearing in S'p then becomes 

[Gp(t,{*-}r' = 2iiG°((;)|-'-E;(k, {ri) 

(fc) + IVF^spiT) 

A relevant quantity to compute is the reduction of the density of states per spin resulting from SP fluctuations 
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FIG. 2. Leading one-particle self-energy correction due to SP fluctuations, for electrons at +fcF. The wiggly (dashed) line 
corresponds to the lattice fleld (—kp) propagator xsp- 




3) 




FIG. 3. Dressed electron-hole bubbles in the a) Landau and b) Peierls (or Cooper) channels, 
(dashed) lines correspond to dressed and bare electron propagators, respectively at +kF {—kp)- 



Thick and thin continous 



(pseudo-gap effect). After analytic continuation of [Gp] ^ to real frequencies, this is found to be 

zi 
ttL 



1 a[2(^ + «]i/2 ^^^^ 



TTVF [2{'d + k) - a'^]'d' 

where the presence of the renormalization factor zi ensures that the pseudo-gap effect is the result of lattice fluctua- 
tions. Following the notation of Rcf. we have 

^=C.((|$-r))"V2 

K = l + -a^ - cj^ 

d={K^+Cj^aY^^ (33) 

where we have defined — z^'^{z~Y I Therefore, as the temperature decreases below the characteristic 

Tgp, lattice fluctuations grow and will progressively freeze electronic degrees of freedom. In turn, this reduction will 
affect the RG flow of various quantities considered in §111. The inclusion of lattice fluctuation effects in the RG flow 
will be considered to lowest order where vertex corrections due to exchange of lattice fluctuations are neglected. For 
example, at one-loop level of the RG, this amounts to substituting G'(fc, {^~}) for one of the propagators appearing in 
the formal expression of the Landau, Peierls, and Cooper elementary susceptibilities (Figure 3 and Appendix A). At 
the two-loop level, all the relevant next-to-leading singular diagrams shown in Figure 4 involve absorption and emission 
of an electron-hole pair at small q in the intermediate state, which can be dressed by lattice fluctuations (Appendix 
B). From the results of Appendix A and B, this amounts to replacing the outer shell logarithmic contribution of all 
diagrams by the following expression 

dl -> TTVF D[EQ{e)/2, {(t)-}]d£, (34) 

where D[Eq{£)/2, {0~}] is density of states ( ^ ) in the presence of a pseudo-gap. 

Substituting this RG generator in the flow equation for the renormalization factor z^^ in (|ll|) will slow down the 
electronic contribution to the decay of the quasi-particle weight below Tgp. Following the example of the purely 
electronic case, the approximate "scaling" (transient-free) form (|l^) then becomes 

z^\T) « z^\TlF){T/TlFf^^^\ (35) 

where zf^T^p) is given by ( pT| ) and 0*(T) ~ nvFD[T, {(f) }]9* becomes a temperature dependent exponent below 
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FIG. 4. Two-loop a) one-particle self-energy and b) four points vertex singular diagrams involved in the RG flow including 
the leading order corrections due to lattice fluctuations. 



3. Staggered response functions 

The pseudo-gap will affect electronic correlations of the (p, M) channel as well. Indeed, substituting (Q) in ([l^ ) 
yields the following approximate "scaling" form 

x;i'(T, {</)-}) «xf(Ts°p)(r/r°p)"^*^^\ (36) 

where xf!{Tsp) is given by (|l^ and 7*(r) = ttvfD[T, {(f)^}]'^*. 

As for the imaginary part (|19|), the above results and Appendix A. 2 lead to 

Imxf (q + 2fcF, LO, {<^-}) = -xf (T, {0-})Imx(q + 2fcf , u, {0"}) 

^_^M(j.{^-})^MWM^^, ,^0 (37) 

which is further reduced by lattice fluctuations. 



^. Magnetic susceptibility 

The calculation of the spin response function Xs at small {q,uj) is known in the purely electronic situation |p^ . Its 
generalization when lattice fluctuations are present is straightforward and yields 

where in the above scheme of approximation ^(q', {4>~}) is the elementary electron-hole bubble dressed by lattice 
fluctuations (Fig. 3a and Appendix A.l), and gi{T) is given by ( ^ at £t- According to the results of Appendix A, 
one finds 

pvpq 



pvpq — ijJ 

p 



Rex(g, a;, {</.-})= / D[u;' , {r}] ( -^j) da;' ■ 
Jo \ 5^ J ^ 

lmx{q,uj,{(f)^})= J D[uj',{(j)^}](^-^^ daj'^pTTVpqSipvFq-uj), (39) 

for the real and imaginary parts of the dressed electron-hole bubble. It follows that in the static (a; — > 0) and uniform 
(g — > 0) limit, the temperature variation of the spin susceptibility Xs{T, {4>^}) will be depressed in the presence of 
lattice fluctuations that grow up below Tgp (Fig. 5. a ). In the very low-temperature domain, the SP correlation length 
becomes exponentially large and the magnetic susceptibility becomes thermally activated. When these temperature 
conditions prevail, the system is almost long-range ordered. 
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FIG. 5. Calculated low temperature variation of a) the uniform magnetic susceptibility and b) nuclear relaxation rate in 
the presence (continuous line) and absence (dashed line) of lattice fluctuations. 

5. Nuclear relaxation rate 

From the above results, the temperature dependence of the nuclear spin-lattice relaxation rate Tj~^ can be readily 
calculated. The basic expression for Tj~^ in one dimension is well known to be 

Tfi =M P T / dQ ^^"^(^'"\ (40) 
J ^ 

which can be taken in the limit w — > 0. Here Im^ is the imaginary part of the retarded spin response function and A is 
a constant proportional to the hyperfine coupling It is well established that uniform {Q ~ 0) and AF {Q 2k p) 
spin fluctuations give the essential contributions to the relaxation in one dimension ||l^ , allowing in turn to make the 
following decomposition 



Tf 1 = Tf i[g - 0] + Tf i[Q - 2kF]. (41) 

For the staggered part T^'^IQ ^ 2fci?], we have seen from ( |37| ) that Imx(2fcF + 9, w) ( whose definition coincides with 
— Imx^^Q(2fci? +(7,a;)) is peaked when q lies in the interval [—T/vp,T/vf], which leads to 

T-\Q ^ 2kF] 2. C,TD[T, {r}]x+(r, {0"}), (42) 

where Ci = nvp^ | A ^ tanh(l/4). 

As for the uniform contribution, the use of (^8[]39|) immediately leads to 



T-\Q ^ 0] ^ CoT ^^^'^^ ^. (43) 

[l~y,iT)D[T,{^~}]]' 



where Cq = 47r(uF)~^ | A p. 

As expected, the reduction of spin fluctuations by short range lattice correlations below Tgp will decrease the 
amplitude of the relaxation via the reduction of the density of states. At suflticiently low temperature, namely when 
the correlation length S,sp becomes exponentially large, one flnds ~ Xs ~ g-pTsp/T ^ ^j-^g high-temperature 
regime, lattice fluctuations are small and the uniform component T^^ CqTxI eventually dominates the relaxation. 
Using the scaling form (^) in ( ^2| ) and (^), the temperature dependence of Tj~^ is summarized in Figure 5.b where 
it is compared to the case without fluctuations. 
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FIG. 6. Comparison between calculated (continuous line) and observed (crosses) temperature profiles of magnetic suscepti- 
bility for (TMTTF)2PF6. The data are taken from Ref.[7], 



IV. APPLICATION TO THE SP SYSTEM (TMTTF)2PF6 



By way of application of the present theory to the sulphur based organic compound (TMTTF)2PFg, Previous anal- 
ysis of this system in the normal phase using the one-dimensional electron gas model can be used for the determination 
of the input bare parameters of the purely electronic part of the model. NMR analysis of Wzietek et al, [|l^ have 
shown that 5 1 ~ 52 — -9, with Ep ~ 1600K, give a rather good description of the temperature dependent magnetic 
susceptibility in the high temperature domain of this material. As for the observed characteristic temperature scale 
Tp « 220K (below which the system presents insulating properties), it can be used to identify, together with (p^), 
the low-temperature domain of strong umklapp scattering, which then allows to take ^3 « .2 As for the input 
parameters for the lattice component, one will fix the value of Tgp at 60K which is the characteristic temperature 
scale for the onset of strong lattice fluctuations in x-ray experiments . From the above set of figures all quantities 
of interest can be calculated. 



A. Magnetic susceptibility 



The temperature-dependent EPR spin susceptibility (TMTTF)2PF6 measured by Creuzet et al, Q is reported in 
the low-temperature domain in Figure 6 . Xs{T) decreases monotonously from the high temperature domain and 
becomes weakly temperature dependent near 80K, which is typical of all members of the sulphur series in the normal 
state. However, in the low temperature domain below 60K, the spin susceptibility decreases by roughly 40 % down 
to the true SP transition at Tsp « 19K, below which it becomes thermally activated. Using the above set of 
parameters for the model, the theoretical prediction for Xs{T, is illustrated in Figure 6 and gives a fairly good 

description of the temperature variation of Xs in the SP pseudo-gap regime. At very low temperature, the present 
theory predicts a thermally activated behaviour when grows exponentially, which is found to mimic the actual 
temperature dependence below the true transition temperature at 19K. However, a more realistic description of the 
SP system in this low temperature region, would require the inclusion of the interchain coupling. 



B. Nuclear Relaxation 



The temperature profile of Tf^ for (TMTTF)2 PFg, measured by Creuzet et al., [§ is given in Figure 7. From the 
analysis of nuclear relaxation in the high temperature domain Tgp < T < Tp, where SP fluctuations are weak, the 
contribution to T-^^ is well known to be purely electronic in character. A quantitative description of the relaxation rate 
in this regime can be obtained from (|4|) and (|43| ) neglecting the dependence on |l^. Below Tgp, the relaxation 
rate shows a 30% decrease between Tgp and Tsp due to one-dimensional lattice fluctuations. According to ( ^2[) and 
( ^) both the staggered and uniform parts of the relaxation are affected below Tgp. Thus, using these expressions 
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FIG. 7. Comparison between calculated (continuous line) and observed (crosses) temperature profiles of nuclear relaxation 
rate for (TMTTF)2PF6. The data are taken from Ref.[7]. 



for the above set of parameters, one obtains the TT^ temperature profile shown in Figure 7. The theoretical curve 
is obtained from the expressions (^3|), (j^) and (36) in which the values of the constants Co{ttvf)~^ ~ 12.1 and 
Ci{TrvF)^^Tpxt{Tp) — 9.6 results from the analysis of T^^ data made by Wzietek at, in the high temperature 



domain T > T^p. 
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APPENDIX A: ELEMENTARY SUSCEPTIBILITIES IN PRESENCE OF PHONONS 

In this appendix, we proceed to the calculation of the elementary electron- hole bubbles dressed by lattice fluctuations 
in the Landau, Peierls and Cooper channels, respectively. 

1. Susceptibility at small q and uj 



The expression corresponding to the electron-hole bubble of Figure 3a is 

T- 



L 



(Al) 



p.k 



for one spin orientation. Using the spectral representation of Gp{k, {(j) }) and performing the fermion frequency sum, 
one finds 



dio'lmGpik,io',{r}) 



n[ep{k + q)] - n[w'] 
icum + £p{k + q) — Uj' 



(A2) 



Since most of the spectral weight appears in the region lo' w ep{k), one can replace ep{k) by lu' in the ratio of the 
above integral, which can be cut off at ±Eq/2. From the definition of the density of states (^2[), one finds, after 
analytic continuation to real frequencies. 
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2j_E„/2 V J ^pvpq-oj 

^ r+Eo/2 _ f 5n\ v-^ 

Imx(g,a;,{<^ }) = - j ^ D[uj',{(j) }] (^-^j (hj'2_,P'^VFq5{pvFq - oj). (A3) 

Here D{uj,{(j)^}) is the density of states per spin in the presence of the spin-Peierls pseudo-gap. The calculation of 
x(9){0~}) at the step £ of the renormalization group procedure gives the same expression, except for which is 
replaced by Eo{£). 



2. Peierls and Cooper susceptibilities 

The calculation of the Peierls electron-hole bubble of Figure 3b starts with the following expression 

{<t>-}) = -rYl G-{k,u;n, {ct>~})Gl{k + Ikp + + oj^), (A4) 

for both spin orientations. At zero external variables and after a summation over the fermion frequencies and the use 
of the spectral representation for G-{k, {<;^~}), one gets 



+°" n[u;'] - n[-€-{k)] 



x(2fcF,{r}) = / cLj'lmG-ik,Lo,{cl>-}) 



uj' + e-{k) 



(A5) 



which actually coincides with the real part of the Peierls bubble. Assuming that ImG'_(fc, w', {0 }) is peaked in the 
region oj' w e-(fc), from which one replaces e-{k) by ui' in the ratio appearing in the integral, we find 

X(2fc^, {</>-}) = - r^'D[u.',{r}W'-^^^^^^. (A6) 
Jo ^ 

As for the electron-electron (Cooper) elementary bubble at zero external variables corresponding to the expression 

Xilr}) = j^^G.{k,{'l>-})Gl{-k), (A7) 

k 

the property G^{—k, —cOn) = —G'^{k + 2A;f, w„) leads to the relation xC^kp, {4'~}) = ~x{{4'~})- 
Within the renormalization group scheme at £, the same expressions become 

dx{2kF, {r }) = -D[Eo{e)/2, {ct>-}]di 

= -dxar}), (A8) 

when evaluated in the outer energy shell. Finally, after analytic continuation to real frequencies of (A4), we can carry 
over the same type of calculation for the imagnary part of the Peierls bubble at T in the limit of small frequency and 
we find 

Imx(2A;F + q,uj^O, {(/)-}) ~ -- / dLj'D[uj', {n[io'] - n[-u;' + VFq])S{ij' + {u - VFq)/2) 

2 J-Eo/2 

8rcosh2(/3t;Fg/4) 



APPENDIX B: ONE-PARTICLE SELF-ENERGY AND FOUR-POINT VERTEX PART IN PRESENCE 

OF PHONONS 

1. One-particle self-energy 

The expression for the one particle self-energy diagram of Figure 4a reads 
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L2 

k',q 

f +00 



where the second Hnc results from a fermion frequency summation and the use of the spectral weight representation. 
The approximation scheme of Appendix A then allows to put e-{k') « w' in the ratio appearing in r.h.s of this last 
expression and to cut off the integral over w' at ±£Jo/2. In the RG procedure, the outer energy shell evaluation of 
S_|_ at £, which is obtained after the frequency sum over Um leads to dS_|_ = dS^ + dS^, where 

= -^D[Eo{i)/2,{r}] dEo{£) 

X /^g l^t^-^offl/^-^fg ] - n[±Eo{£)/2]){nB[-VFq]+n[vFq + £+(&)]) 



Sttvf 



2vFq + ioJn - e+(fc) 

[G^(k)]-^D[Eo{t)l2,{r}]dt, (B2) 



to leading order in [G"(fc)]^^ = zcj„ — e+(fc), for /3Eq{£) ^ 1, and where ^^[.x] = (e^^ — Here the integration 

over the momentum transfer q is found to contribute only in the interval 2vFqo > 2vf \ q \ > Eo{£), where qo is a 
momentum transfer cut-off. 



2. Vertex part 



The evaluation of the two-loop vertex part in the presence of lattice fluctuations proceeds along similar lines. The 
diagrams of Figure 4b corresponds to ri^^ + T^^ and to the generic expression 

^aiar}) = -9'^J2 r )G° (fc' - ?)G^(fci,2 + ?)Gi(^i,2 + q'-q). (B3) 

In the RG sense one can drop the dependence on the external variables {ki, k^, q}. Performing frequency sums on 
and LUm' and using the above approximation scheme, one gets 

d^al ^ 9'^^D[Eo{£)/2, {</,-}] dEo{£) 

X J ^(^[n[Eo{£)/2 + VFq']-n[Eo{£)/2] [nB[-VFq'] + n[vFq']] + Eo{£) ^ -Eo{£)) 

for 2vFqo > 2vf \ q' \ > Eq{£), (3vFq' » 1, and PEq{£) > 1. 
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